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FOREWORD 

This paper was presented, upon 
invitation of the symposium chairman, 
as part of the Beckmann Award Sym- 
posium at the 131 st meeting of the 
American Chemical Society on £ April 



1957 in Miami, Florida 



ABSTRACT 



Programming of the operational 
analogue computer for the rapid 
solution of problems dealing with 
complex chemical reaction-rate 
mechanisms is described. As an 
illustration, an enzyme-catalyzed 
reaction is considered and it is 
shown that computer results agree 
well with experiment, but not with 
the steady-state approximation 
customarily made for this system. 

It is pointed out that the oper- 
ational analogue computer is a 
relatively simple and inexpensive 
set-up, and that programming prin- 
ciples, as developed and described 
here, make it applicable to a wide 
variety of theoretical and industrial 
problems. 



. 






An operational analogue computer is an assemblage of 
independent high-gain DC amplifiers, with accessory resistors, 
capacitors, and potentiometers, and a peg-board arrangement 
for making connections. Each time-dependent variable is 
represented on the computer as a voltage, and the amplifiers 
can either add together several such voltages, or else add 
and integrate. It is the purpose of this paper to show how 
an analogue computer of this type can be used for investigating 
the mechanism of a complex chemical reaction. The unravelling 
of such a mechanism would then, for example, permit ready 
scale-up from laboratory to plant size equipment, and thus be 
of valuable assistance in the selection of optimum size reactors. 

In many actual applications of chemical kinetics, the rate 
law corresponding to a given mechanism is approximated by the 
use of simplifications suitable to the problem at hand - such 

as the steady-state hypothesis customarily required in a large 

3 - • 

number of chemical mechanism studies. When, however, the 

concentration.of each species (including intermediates) is 
desired as a function of time, especially in the early stages 
of reaction before the steady state is established, it becomes 
necessary to use the complete integrated form of the rate law. 
This form is also required if the rates of several steps are 
of comparable magnitude, so that the concept of a rate- 
determining step is not significant. These complete expressions 
are extraordinarily complex and, indeed, in many cases can be 
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evaluated only by numerical methods. This is especially true 
for reactions involving bimolecular or trimolecular steps in 
the mechanism. 



An alternative treatment for the complete rate lav/ is that 
of working with the set of simultaneous differential equations 



directly obtainable from the mechanism. Each of these 
equations is of first order and first degree, though non- 
linear for bimolecular or trimolecular steps. (It should be 
noted that even though the reaction is of higher than first 
order in the chemical sense, each differential equation is 
mathematically of the first order.) 

Consider, for example, a reaction of two consecutive 
unimolecular steps; 



Selecting the rate equation for component B as an example, 
and using lower-case letters to represent concentration, 



6T 




( 1 ) 



db 

dt 



ki a - k^b 



( 2 ) 



Defining the operator 1/P by the equation 




0 



( 3 ) 



equation 2 may be integrated to 



t 




(4) 



2 - 






i ■ 

- 

• I 



3fV M »•#-! 








' 






FIG. I THE INTERGRATING AMPLIFIER! (upper) 
ELECTRICAL CIRCUIT; (lower) SCHEMATIC. 
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or 



b 



( 5 ) 



- -1/P (k 2 b - k^a) + b Q 

in which subscript zero refers to the concentration at 
zero time. 

The upper portion of Fig. 1 shows the electrical circuit 
for an amplifier arranged to perform the operation indicated 
by equation 5; this is shown schematically in the lower 

7 

portion of the figure, using conventional computer symbols. 
The amplifier, with amplification -A, is shown as a triangle, 
the grid input being to the base of the triangle, the plate 
output from the apex. The potentiometers, set at a^and a 2 , 
control the fraction of the input voltages, e^and e 2 (referred 
to ground) fed through the input resistors, R^and R 2 , to the 
grid of the amplifier. For feed-back capacitance C and with 
output voltage e, (also relative to ground) it follows from 
Kirchhoff’s Law, applied to the grid input, that 



a e - e 
11 g 

R, 



a e - 
2 2 



JS + C 



R, 



d(e - e) 

s. ♦ i g - 0. 

dt 



( 6 ) 



where subscript g refers to the grid. Now the grid voltage 
and current are very small for a well-designed high-gain 

5 7 

amplifier, and from this consideration it develops* * that, 
in terms of the operator l/P, the output voltage becomes 



e = -1/P( 



TTTT 



R 2 C 



2 ) 



+ E, 



( 7 ) 



where E Q is the initial value of the output voltage 
(disconnected when computer operation begins). 
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FIG. 2 THE SUMMING AMPLIFIER (upper) 

SET-UP FOR THE ENZYME PROBLEM, (lower) 
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The similarity in form between equations 5 and 7 suggests 



the use of such amplifiers in solving the type of problem 
exemplified by equation 1, The requirement is that the vol- 
tages suitably represent concentrations and the resistor- 

capacitor-potentiometer combinations represent the rate 
♦ 

constantSo 

In many problems it is also necessary to add together 
variables. For this purpose can be used the summing am- 
plifier, as shown schematically in the upper portion of 
Figure 2, The arrangement is similar to that for the inte- 
grating amplifier, save that the feed-back is resistive, of 
value R^. Making assumptions as previously, the equation 
for the output voltage is 



For single input and equal resistors, the summing amplifier 
functions as a sign changer. 

If the mechanism under consideration contains steps of 
molecularity greater than one, it is also necessary to use 
a function multiplier, which is an electronic device to 
generate the product of two variable voltages. For commer- 
cially available multipliers the output, as a function of 
input voltages e^and & 2 > is given by 



a 




e 



( 5 ) 



e 



- me,e 



12 ’ 



(9) 
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where m is a constant of proportionality (equal to 0.01 for 
the equipment used in this laboratory). On the schematic 
diagrams the multiplier is shown as a rectangle with an X 
at its center. Another auxiliary piece of equipment is a 
function generator, whose output conforms to some arbitrarily 
given value for any given input. This device is useful for 
non-isothermal reactions, where rate constants may change. 

To illustrate the use of computer technique, the following 
mechanism has been selected: 

k k 

A + B 1 w C 3 ^ B + D (10) 



in which B is a catalyst of some sort. Equation 10 could, 

for example, represent the Michaelis-Menten formulation for 

4 

enzyme activity, B symbolizing the enzyme and A the sub- 
strate. Assuming only A and B present initially, the 
operator equations for A, C, and D become 



a = - 1/P 

-c = - 1/P 

d = - 1/P 



k n ab + k^(-c) + a 



] 



0 



k^b + k 2 (-c) + k^(-c) 



k^(-c) 



J 



(ILL 

( 12 ) 

(13) 



The equations are written in terms of -c to avoid the use 
of additional amplifiers for sign changing. The equation 
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for B can be deduced from the stoichiometry of the system; 



-b 



0 




(14) 



-b being used here since the function multiplier changes sign. 
These four equations suggest the use of at least four 
amplifiers (three integrating and one summing) plus a function 
multiplier. The arrangement used to solve this problem is 
shown schematically in the lower portion of Figure 2. 

The specific problem chosen for illustration was the 
enzyme-catalyzed oxidation of an appropriate oxygen acceptor 
by hydrogen peroxide, a system which has been studied by 
Chance . Typical values of the parameters for this problem 
are: 

a -6 -6 7 

0 = 4 x 10 M, b Q = 1x10 M^k-,^ = 1x10 , k 2 = 0.2, 1^=5 



(t in seconds for the rate constants). To illustrate the scaling, 
it will suffice to consider only the amplifier generating 
a, the voltage representing a. Let = a/a,<Xg - b/b 

£ = c/c, where each V is a magnitude-scaling factor, and 

a", ^T, and c the voltages representing the respective con- 
centrations. Then equation 11 may be rewritten 



a = 



1/P 



k l ^ B &b 



+ k' 



<^C 

<^A 



(-c) 



0 



(15) 



To take care of time-scaling, a factor is introduced with 
each feed-back capacitance, such that ft seconds on the 
computer corresponds to one second of actual time. Then by 



applying the electrical equations 7 and 9 to this amplifier, 



a - -1/P 

where subscript 1 refers to the portion of the circuit into 
which mail) is fed, subscript 2 to that for -c. By equating 
coefficients in equations 15 and 16, 






ab 



R-jC 



+ ^2 

r 2 c 



£ 



(- 0 ) 



+ a 



0 



( 16 ) 



k 



1 




RjC 



k 2 = a 2 P 

r 2 c 



( 17 ) 

(IS) 



These two equations, combined with similar equations for 
the other amplifiers, permit the proper selection of the 
resistors, capacitors, and potentiometer settings. In Table 
1 are shown values of these components corresponding to one 
case considered, and in Figure 3 are shown simultaneous 
traces (from a four-channel recorder) of each of the concen- 
trations vs time for this case. Concordance with the results 
of Chance are quite satisfactory; e_,£, , the maximum 
concentration of the intermediate C and the time required 
to attain this maximum agree quite well with the results 
reported by him. 

The assumption that there is established a steady state 

4 , 

in b and c, after an initial transient period leads to 
the relationship 
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Table I. Typical Values of Components for the Enzyme Problem 



Amplifier (Species) 





A 


B 


c 


D 


Magnitude-scaling factor, 
cK , x 10 d 


10 


1 


1 


10 


Feed-back 


1/f 


.2 meg 




v< 


Input Resistors, megohms 


1, 10 


.2 o 1 , 


1, 1 


1 


Initial voltage 


3S 


95 


0 


0 


Initial concentration, 
micromoles/liter 


3.S 


«95 


0 


0 



Potentiometer (Rate-constant) 

3 

0o5 
5 



Time-scaling f actor, p , =10 

Multiplier scaling factor, m, = 0„01 



Potentiometer Setting 
Rate Constant (t in sec®) 



1 

1.0 

1 x 10' 



2 

0.02 

0.2 
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CONCENTRATION IN MICROMOLE PER LITER 




t,sec. 



FIG.3 



KINETICS OF THE ENZYME PROBLEM. 




FIG. 4 EFFECT OF CATALYST CONCENTRATION 
ON STEADY— STATE ASSUMPTION. 



12 



c 



(19) 



k 

2 



k 

+ 3 



k 



1 



ab, 



which applies strictly only when db/dt = dc/dt = 0, In 
Figure 4 is shown a portion of the plot of c v£. ab, for a^ 

= 4 x 10 ^ m, and various values of b^, as made directly 
from the computer on an X-Y recorder, (note that time 
increases from right to left in this figure.) The straight 
line is that computed from equation 19» It can be seen that 
if bp is comparable with a^, the steady-state approximation 
is not valid, even near the completion of reaction. Only at 
small values of b^/a^ do the actual curves approach equation 
19, although the straight line could serve as a rough average 
for the latter stages of reaction. Whereas deductions such 
as these are readily made from computer traces, it is indeed 
quite difficult, if. not impossible, to produce them by 
ordinary mathematical methods® 

For a general kinetics problem at least one amplifier 
will be required for each component participating in the 
mechanism, plus such others as may be required for sign- 
changing; one potentiometer for each rate-constant; and one 
function multiplier for each biraolecular step (or two for a 
trimolecular step, etc.). The speed, simplicity, and 
adaptability of the operational analogue computer makes it 
also a useful device for the determination, from observed 
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kinetics data, of the rate constants for assumed mechanisms. 
Trial values of the various rate constants are used, by 
changing potentiometer settings and, if necessary, resistors, 
until a fit with the data is obtained. Since changing 
parameters ordinarily takes only a few seconds, such a 
problem can be solved quite rapidly. 

It follows that the independent variable of the problem 
need not be time, but could also be any quantity varying 
linearly with time, using appropriate scaling. It is thus 

possible to study the kinetics of steady-flow processes. 

2 

An example is the study of catalytic isomerization ; the 
results of mechanisms assuming either single-site or dual 
site adsorption are readily obtained for a wide variety of 
initial conditions and flow rates. 

A computer of the sort described is of only moderate 
expense - a quite satisfactory combination need cost no more 
than a good polarimeter - and can be operated by a single 
chemist or chemical engineer working in his own laboratory. 

A problem of the type described can be set up and solved in 
a few hours. 
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